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Abstract. We present a self-consistent, absolute isochronal age scale for young (< 200 Myr), 
nearby (< 100 pc) moving groups, which is consistent with recent lithium depletion boundary 
ages for both the /? Pic and Tucana-Horologium moving groups. This age scale was derived using 
a set of semi-empirical pre-main-sequence model isochrones that incorporate an empirical colour- 
Teft relation and bolometric corrections based on the observed colours of Pleiades members, with 
theoretical corrections for the dependence on logp. Absolute ages for young, nearby groups are 
vital as these regions play a crucial role in our understanding of the early evolution of low- 
and intermediate-mass stars, as well as providing ideal targets for direct imaging and other 
measurements of dusty debris discs, substellar objects and, of course, extrasolar planets. 
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1. Introduction 

Over the past couple of decades several hundred low- and intermediate-mass stars have 
been identified within ~ 100 pc of the Sun as part of a concerted effort to identify and 
characterise the young solar neighbourhood. The distribution of these stars on the sky 
is not uniform, instead they comprise dispersed, predominantly unbound associations 
in which the members share a common space motion. Although the age ordering of 
these young associations is well-constrained i.e. the AB Dor moving group is older than 
the Tucana-Horologium moving group (Tuc-Hor), which in turn is older than the /3 Pic 
moving group (BPMG), the absolute ages of these groups are still under-constrained. 

Theoretical model isochrones are (arguably) the most commonly used method of age¬ 
dating young (presumably coeval) stellar populations. Although the use of maximum- 
likelihood fitting techniques to subjectively fit model isochrones has been used for more 
distant Galactic clusters (e.g. NGC 2547 in Naylor & Jeffries 2006), such methods have 
not yet been applied to the young groups within ~ 100 pc. Additionally, for these young 
groups there has been a distinct lack of homogeneity when it comes to the fitting of model 
isochrones, particularly in terms of the models adopted and the photometric bandpasses 
used to construct the colour-magnitude diagrams (CMDs). Furthermore, recent lithium 
depletion boundary (LDB) ages, which are advocated in the recent review of young stellar 
ages by Soderblom et al. (2014) to provide our best chance of establishing a reliable and 
robust age scale in the range 20 — 200 Myr, have been demonstrated to be systematically 
older when compared to isochronal ages (see e.g. Kraus et al. 2014). 

In this contribution we provide a self-consistent, absolute isochronal age scale for young 
(< 200 Myr) moving groups within 100 pc of the Sun. This age scale is based on homo¬ 
geneous fitting of photometric data in the My, V — J CMD using a maximum-likelihood 
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fitting statistic in conjunction with semi-empirical model isochrones. We find that the 
isochronal ages for both the BPMG and Tuc-Hor are consistent with recent LDB ages, 
thereby instilling confidence in our isochronal age scale. 

2. Sample of young, nearby moving groups 

In this contribution we focus on the young, nearby moving groups within 100 pc; specif¬ 
ically the AB Dor moving group, Argus association, BPMG, Carina association, Columba 
association, r) Cha cluster, Tuc-Hor, TW Hya association (TWA), and 32 Ori group (see 
also the recent review by Torres et al. 2008). For each group, we have assembled a list 
of members and candidate members (hereafter simply referred to together as ‘members’) 
from the literature. For the inclusion of candidate members, we require that they have a 
membership probability of ^ 90% as calculated using the Bayesian Analysis for Nearby 
Young AssociatioNs (BANYAN; see e.g. Malo et al. 2013). Furthermore, we require that 
any candidate members have a measured radial velocity and/or trigonometric parallax 
which are/is consistent with membership for a given young group. 

Our sample of members (not accounting for unresolved multiples) includes 89 members 
of the AB Dor moving group, 27 members of Argus, 97 members of the BPMG, 12 
members of Carina, 50 members of Columba, 18 members of 77 Cha, 189 members of 
Tuc-Hor, 30 members of TWA, and 14 members of 32 Ori. 

Whilst all of our members have counterparts in the Two-Micron All-Sky Survey Point 
Source Catalog (2MASS PSC; Cutri et al. 2003), near-IR CMDs are not ideal for age 
determination primarily because for 2500 < T e ff < 4000 K the loci of young clusters 
becomes vertical (J — K s ~ 0.9mag), and hence degenerate with age. Therefore we 
supplement our near-IR photometry with P-band data and derive ages for our sample 
of groups in the My, V — J CMD so as to minimise the effects of circumstellar material 
on the K s magnitudes for the youngest groups (especially 77 Cha and TWA). 

For assigning distances to each star in our list of members we prefer to adopt trigono¬ 
metric parallax measurements, however when this is not an option (or where the parallax 
estimates appear to be in error or have unusually large uncertainties) we either adopt a 
kinematic distance from the literature or derive one using the ‘moving cluster’ method 
(see e.g. Mamajek 2005). 

3. Semi-empirical model isochrones 

There are large discrepancies between pre-main-sequence (pre-MS) tracks and obser¬ 
vational data (such as dynamical masses, age trends, etc.; see e.g. Hillenbrand & White 
2004; Feiden this volume). In Bell et al. (2012) we discussed the reasons why the Pleiades 
represents an ideal benchmark cluster for the comparison of the models to the data, how¬ 
ever chief among these is the fact that both the distance and age have been determined 
independently of the use of model isochrones. Fig. [TJ which shows the V,V — I c CMD of 
the Pleiades with several commonly adopted sets of model isochrones overlaid, effectively 
demonstrates the aforementioned discrepancies. At T e g < 4300 K, the models overesti¬ 
mate the flux in the optical by up to a factor of two. Hence, if we are to derive consistent 
ages for young stellar populations using pre-MS model isochrones we must first perform 
some form of empirical correction to the models for T e s lower than ~ 4300 K. 

In Bell et al. (2013, 2014) we introduced a method of creating semi-empirical pre-MS 
model isochrones using the observed colours of young stars in the Pleiades, in addition 
to incorporating theoretical corrections for the dependence on the surface gravity (logg). 
Briefly, we used binary systems with well-constrained dynamical masses to demonstrate 
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Figure 1 . The V, V — I c CMD of the Pleiades with several sets of model isochrones overlaid. 
The upper reference refers to the interior models, whereas the lower reference corresponds to 
the colour- T e s relation and bolometric corrections adopted to transform the model into CMD 
space. 

that the model predictions for the A' s -band are essentially correct. This then allows 
us to use the AVband magnitude along our fiducial locus as a T e g indicator, but also 
ensures that our semi-empirical model isochrones have a reliable mass scale tied to that 
of the fiducial cluster. On the assumption that the A s -band magnitudes predicted by the 
models are correct, we can then use CMDs with the A' s -band on both axes (e.g. K s , V — 
Ag) to calculate any discrepancy between the fiducial locus and the model isochrone, 
which in this case would correspond to the requisite empirical correction in the P-band 
necessary to fit the Pleiades at a given age and distance. This empirical correction is 
then applied to the theoretical bolometric correction grid (calculated from atmospheric 
models and a function of T e g and logg) at the appropriate T e g irrespective of its logg. 
This process is then repeated for other photometric bandpasses to create sets of semi- 
empirical bolometric corrections. The semi-empirical pre-MS model isochrones used in 
this contribution are based on existing stellar interior models coupled with these newly 
derived semi-empirical bolometric corrections. In this contribution we adopt four sets of 
interior models, namely those of Dotter et al. (2008), Tognelli et al. (2011), Bressan et 
al. (2012), and Baraffe et al. (2015; hereafter referred to as Dartmouth, Pisa, PARSEC, 
and BHAC15 respectively]!] 

4. Isochronal age scale for young, nearby groups 

We use the r 2 fitting statistic of Naylor & Jeffries (2006) and Naylor (2009) to derive 
ages from the My, V — J CMDs of our sample of young groups. The r 2 fitting statistic 

f A subset of these models is available via the Cluster Collaboration isochrone server 
http://www.astro.ex.ac.uk/people/timn/isochrones/ 
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Table 1 . Ages for the young groups in our sample. The penultimate row lists our final adopted 
age for each group for which the associated uncertainties represent the statistical and systematic 
uncertainties added in quadrature.The final row shows literature LDB ages for the BPMG and 
Tuc-Hor (see Section p] for references) which highlights the consistency between the two age 
diagnostics. 


Model 




Group age (Myr) 





AB Dor 

Argus 1 

BPMG 

Carina 

Columba 

77 Cha 

Tuc-Hor 

TWA 

32 Ori 

BHAC15 

145+889 

69±f 

25 ± 1 

46±f 

44±® 

12 ± 1 

46^ 

10 ± 1 

23+2 

Dartmouth 

135+“ 

60±« 4 

23 ± 1 

«±5 3 

43±« 

10 ± 1 

48^ 

7-1 

20±i 

PARSEC 

151 tS 

6Q + 268 

25 1\ 

45±( 2 

43±| 

14 ± 1 

46 ± 1 

13 ± 1 

25t 2 5 

Pisa 

166 tit 

55±f 5 

20 ±1 

41+ 8 

38 ± 3 

8± 1 

40^ 

9 ± 1 

18 ± 1 

Adopted 


- 

24±3 

45+“ 

42+® 

11±3 

45±4 

10±3 

22 1% 

LDB age 

- 

- 

24 ±5 

- 

- 

- 

40 ± 3 

- 

- 


Notes: 

x No final adopted age is given as it is unclear whether the stars in our list of members represent a single 
population of coeval stars (see Section |5j). Note also that the PARSEC models only provide an upper limit on 
the age of the association and therefore we do not provide a lower age uncertainty in the table. 


is particularly well-suited to this task because it allows for the effects of binarity, yields 
reliable uncertainties on the derived parameters and provides a goodness-of-fit test. We 
have made some modifications to the statistic since those of Naylor (2009), however for 
the purposes of this contribution we refrain from detailing these and instead refer the 
reader to Bell et al. (2015). 

In essence the r 2 fitting statistic can be viewed as a generalisation of the y 2 statistic 
to two dimensions on the basis that both the model isochrone and photometric data are 
two-dimensional distributions; specifically uncertainties in colour and magnitude for the 
data and the widening of the model isochrone due to the effects of binarity. As in Naylor 
& Jeffries (2006), we define the r 2 statistic as 


t 2 = — 2 In / / Ui(c — Ci,m — rrii) p(c,m) dedm, (4-1) 

i=l, n ■' 

where [/,; represents the two-dimensional uncertainty function and p the probability dis¬ 
tribution of the expected model in CMD space. We create our model distributions using 
a Monte Carlo method to simulate 10 6 stars over a given mass range and populate these 
using a broken power law mass function. For the inclusion of binaries we assume a uni¬ 
form fraction of 50%. Note that the best-fit age is remarkably insensitive to the adopted 
binary fraction. Our full grid of models range from log(age) = 6.0 — 10.0 in steps of 
Alog(age) = 0.01 dex. For each model/data combination, we then calculate the resultant 
t 2 until we find the model which minimises the above function. This is equivalent to 
maximising the ‘overlap’ between the data points and the model. Fig. [2] illustrates an 
example of the best-fitting My, V — J CMDs of Tuc-Hor, and Table [l] shows the best-fit 
ages for each young group in our sample. 

As part of our modification to the statistic, we now assign prior membership proba¬ 
bilities to each individual star in a given catalogue. Consequently, our best-fitting model 
returns posterior membership probabilities. For a given group in our sample, we expect 
a negligible age spread (or equivalently luminosity spread), and hence we can then use 
these posteriors to identify stars which appear to be non-members (based solely on CMD 
position) in an effort to further refine the membership lists of these young group. 
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Figure 2. Best-fitting My,V — J CMDs of Tuc-Hor. The red circles represent fitted data, 
whereas the blue squares denote objects which are removed prior to fitting as they lie outside the 
area of CMD space covered by the grid of models. Top left: BHAC15. Top right: Dartmouth. 

Bottom left: PARSEC. Bottom right: Pisa. 


5. Discussion 

To assess the reliability of our self-consistent, absolute isochronal age scale, we must 
first compare our derived ages to what are considered well-constrained ages for the same 
groups. Model-independent methods (such as kinematic ‘traceback’ or ‘expansion’ ages) 
should offer such a comparison, however studies have demonstrated that the ages inferred 
from kinematic information alone are simply too unreliable and even unreproducible when 
similar analyses have been performed using improved astrometric data (see e.g. Mamajek 
& Bell 2014). LDB ages have been argued to be both accurate and precise to just a few 
Myr (see the discussion in Soderblom et al. 2014), and have recently been calculated for 
both the BPMG and Tuc-Hoi[|] The LDB in the BPMG has been identified by both Binks 

j Recent analyses (see e.g. Somers this volume) have demonstrated that LDB ages may be 
underestimated by ~ 10 — 20% as a result of neglecting the effects of starspots in the evolution 
of pre-MS stars. 
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& Jeffries (2014) and Malo et al. (2014b), who derive an age consistent with 24 ± 5Myr. 
Note that our age for the BPMG is in excellent agreement with the 22 ± 3 Myr derived 
by Mamajek & Bell (2014) which was based on an analysis of the A-, F- and G-type 
members of the group. In addition, Kraus et al. (2014) identified the LDB in Tuc-Hor 
and derived an age commensurate with 40 ± 3 Myr. Our isochronal ages for both of these 
groups are consistent with the LDB ages. 

Of the remaining groups in our sample, our isochronal ages for the AB Dor moving 
group, Carina, Columba, rj Cha, and TWA are all consistent with previous estimates 
in the literature, however the advantage of our scale is the self-consistency which was 
previously lacking. The 32 Ori group has only recently begun to be investigated (see 
e.g. Mamajek 2007; Shvonski et al. 2010) and hence our isochronal age represents the 
first definitive age for this group. Finally, the reader will note that we do not provide a 
final adopted age for Argus, and there are two primary reasons for this. First, if Argus 
is indeed ~ 40 Myr-old, then the group of 5 A-type stars (all of which were proposed as 
members by Zuckerman et al. 2011 and none of which are unresolved binaries) should be 
located on the zero-age main-sequence (ZAMS). The fact that none of these stars are (4 
are over-luminous and 1 is under-luminous with respect to the ZAMS) suggests that these 
stars are not coeval, but instead represent stars at different evolutionary stages. Second, 
of the K- and M-type stars, less than two thirds have CMD positions commensurate with 
an age of ~ 40 Myr, whereas the remainder appear significantly fainter and thus older. 
We are therefore of the opinion that either our list of members for Argus suffers from a 
high level of contamination or that the association is not physical. 


6. Summary 

We present a self-consistent, absolute age scale for eight young (< 200Myr), nearby 
(< 100 pc) moving groups in the solar neighbourhood. Further details concerning the 
methodology presented in this contribution can be found in Bell et al. (2015), but in brief 
our analysis involves homogeneously fitting the My, V — J CMDs to derive best-fit ages 
using semi-empirical pre-MS models isochrones in conjunction with the r 2 maximum- 
likelihood fitting statistic. Our isochronal ages are consistent with currently available 
LDB ages. Such consistency instills confidence that our new isochronal age scale for 
young, nearby groups is robust and hence we suggest that these ages be adopted for 
future studies of these groups. At the moment, we are uncomfortable assigning a final, 
unambiguous age to Argus as it remains unclear whether our membership list for this 
group constitutes a single population of coeval stars. 
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Discussion 

G. Chabrier: More of a (nasty) comment than a question, sorry. I’m not sure I believe 
your age scale given the methodology of creating the semi-empirical model isochrones you 
have described. First, some of the models you have adopted in your analysis have internal 
inconsistencies with regard to the interior physics (e.g. the stellar interior/atmospheric 
boundary condition) and should therefore not be adopted in the first place. Second, you 
derive semi-empirical T e ff-bolometric correction relations at the age of the Pleiades and 
then apply these at younger ages, however there is no reason to believe that such a 
relation will be valid in this age regime. Third, your models rest almost entirely on a 
well-constrained age (and distance) for the Pleiades. 

C. Bell: I appreciate your concerns, however how are we to determine the ages of young 
stellar populations if none of the models fit the data? Given this discrepancy the only 
option is to make some sort of empirical correction to the models. I’ll answer your points 
in reverse order. First, the Pleiades does have a well-constrained age of ~ 130 Myr which 
has been derived from both the high- and low-mass population using distinct techniques 
and methods that rely on different aspects of stellar evolution, and so I believe the age 
we adopt is robust. Second, true we do derive the T e fj-bolometric correction relation at 
130 Myr, however we attempt to account for the log g dependence using information from 
the atmospheric models; specifically how the bolometric corrections vary as a function of 
logg. Third, again true, the majority of models do suffer from internal inconsistencies, 
however even if we adopt those with fully consistent treatments (such as the most recent 
BHAC15 models), there is still a large discrepancy between the models and the data for 
the Pleiades at low T e g. Furthermore, any internal inconsistencies are effectively ‘tuned 
out’ as a result of the empirical corrections we derive and apply to the theoretical T e g- 
bolometric correction relation, and so these should not have a significant effect on the 
derived ages. 

B. ZUCKERMAN: I’m surprised at the older ages you derive for Argus. I respect Beto 
Torres and collaborators enough to find their age of ~ 40 Myr hard to reconcile with 
those you derive. 

C. Bell: I agree, it is perturbing. Interestingly, if we take the Argus members from 
Torres et al. (2008) and run these through the BANYAN analysis tool we find that 



C. P. M. Bell et al. 


only 11 of the 29 stars (38%) appear to be high-probability 90%) candidate members. 
Furthermore, of these, several have calculated kinematic distances which differ from those 
given in Torres et al. by more than 2 — 3er. Whether this is highlighting an underlying 
issue with the Torres et al. membership list or the BANYAN kinematic distances remains 
unclear. 

J. Gagne: On a related note, are the members you adopt for Argus taken directly from 
the BANYAN analyses of Malo et al.? 

C. Bell: Yes, our membership list for Argus is comprised from the studies of Malo et al. 
(2013, 2014a,b). However, we only include candidate members if they have a measured 
diagnostic (e.g. radial velocity or trigonometric parallax) which is consistent with Argus 
membership and do not include the other high-probability candidate members with no 
such measured diagnostic. 



